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ABSTRACT 

Context. The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence 
gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to 
be important in the solar and stellar dynamo processes. 

Aims. In this work, we have explored the effects of the turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model 
with a solar-like rotation law. 

Methods. In a first step, only vertical pumping has been considered through the inclusion of a radial diamagnetic term in the induction 
equation. In a second step, a latitudinal pumping term has been also included and then, in a third step a near-surface shear has been 
switched on in the model. 

Results. The results reveal the importance of the pumping mechanism for solving current limitations in mean field dynamo modeling 
such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. In the case that a meridional flow is assumed to 
be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle 
and the latitudinal distribution of the sunspots activity. In models that consider shear near the surface, a second shell of toroidal field 
is generated above r=0.95R o at all latitudes. If the full pumping is also switched on, the polar toroidal fields are efficiently advected 
inwards, and the toroidal magnetic activity survives only at the observed latitudes near the equator. With regard to the parity of the 
magnetic field, only models that combine turbulent pumping with near-surface shear always converge to the dipolar parity. 
Conclusions. This result suggests that, under the Babcock-Leighton approach, the equartorward motion of the observed magnetic 
activity is governed by the latitudinal pumping of the toroidal magnetic field rather than by a large scale coherent meridional flow. 
Our results support the idea that the parity problem is related to the quadrupolar imprint of the meridional flow on the poloidal 
component of the magnetic field and the turbulent pumping positively contributes to wash out this imprint. 

Key words. Sun: magnetic fields - Sun: activity 



1. Introduction 

Flux-dominated Babcock-Leighton (FDBL) solar dynamo are 
mean field models where the poloidal field is generated at the 
surface by the transport and decay of bipolar magnetic regions 
(BMRs) which are formed by twisted buoyant magnetic flux 
ropes. For this process to occur, differential rotation must be 
able to develop intense toroidal magnetic fields either at the 
tachocline or at the convection zone. Numerical simulations have 
shown that magnetic flux tubes with intensity around 10 4 -10 5 G 
are able to become buoyantly unstable and to emerge to the sur- 
face to form a bipolar magnetic region with the appropriate tilt, 
in agreement with the Joy's law. One important limitation of the 
scenario above is that 10 5 G results an energy density that is an 
order of magnitude larger than the equipartition value, so that a 
stable layer is required to store and amplify the magnetic fields. 
This brings another problem with regard to the way by which the 
magnetic flux is dragged down to deeper layers. 
In the lack of accurate observations of the flow at the deeper lay- 
ers, numerical simulations have shown that the penetration of the 
plasm a is restricted to only a few kilometers below th e overshoot 
layer dGilman & Mieschlr2004t iRudiger et alj|2005l) . neverthe- 
less, the magnetic fields can be transported, not only downwards, 
but also longitudinally and latitudinally when strong density 



and turbulence gra dients are present in the medium due to the 
turbu lent pumping dZiegler & Riidigerf2003HDorch & Nordlundl 
l200!h . 



In axisymmetric mean field models of the solar cycle, the effects 
of the turbulent pumping have been rarely considered. A first ap- 
proach showin g the importance of the pu mping in the solar cycle 
was made by iBrandenburg et all dl992l) . since then few works 
have incorporated the diamagnetic pumping component in the 
dynamo equation as an extra diffusive t erm which provides a 
down ward velocity dKiiker et alj 1200 U iBonanno et all 120021 
120061) . More recently, Kapvla et al. d2006bb have implemented 
simulations of mean field dynamo in the distributed regime, 
including all the dynamo coefficie nts, previously evaluated 
in magneto-convect ion simulations dOssendrijver et all 120021 : 
iKapyla et alj l2006al) . They have produced butterfly diagrams 
that approximately resemble the observations. However, to 
our knowledge no special efforts have been done to study 
the pumping effects in the meridional plane (i.e., inside the 
convection zone) or in a FDBL description. The latter has been 
found to be particularly successful at reproducing most of the 
large s cale features o f the s ol ar cycle (Dikpati & Charbo nneaul 
1 19991: iDikpati et al.1 12004 iGuerrero & de Gouveia Dal Pinol 
|2007a|E hereafter GDPa, b). 
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Fig. 1. Radial and latitudinal profiles for a (continuous line), r\j 
(dotted line) and for the pumping terms y r and yg (dashed and 
dot-dashed lines, respectively). All the profiles are normalized 
to their maximum value. 

In this work, we explore the effects of the turbulent pumping 
on a FDBL model. In a first approximation, we include the ra- 
dial turbulent diamag netism velocity term in the indu ction equa- 
tion as described by iKichatinov & Ruedigerl (1199 23, and then 
in a second approach we add the pu mping terms calculated in 
local magneto-conve ction simulations dOssendrijver et a02002t 
iKapvla et al.l l2006al) . This latter approximation includes not 
only the radial but also the latitudinal contribution of the 
pumping. Finally we will also discuss the implications of the 
pumping when the ne a r-surfa ce radial shear layer reported by 
ICorbard & Thompsonl (2002) is considered. In the next, we 
briefly present the model, our results and then, we outline the 
main conclusions of this work. 

2. The model 

Our model solves the mean field induction equation: 

— = Vx[UxB + £~77 r VxB], (1) 
ot 

where U=u p + fir sin 9 is the observed velocity field, fi is the an- 
gular velocity, B=V x (Ae^) + B^e^ are the poloidal and toroidal 
components of the magnetic field, respectively, tj t is the mag- 
netic diffusivity and 

£ = aB + y X B, (2) 

corresponds to the first order terms of the expansion of the elec- 
tromotive force, u x b, and represents the action of the small- 
scale fluctuations over the large scales. The coefficients of (|2|i 
are the so called dynamo coefficients. Normally, the mean field 
models do not consider the second term on the right hand side 
of the equation, i.e, the turbulent pumping. The first term cor- 
responds to the alpha effect that has been considered in several 
ways in the literature. A pure Babcock-Leighton model is an aQ. 
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Fig. 2. (a) Butterfly diagram and latitudinal snapshots for the 
toroidal (b) and the poloidal (c) fields. The dark (blue) and light 
(red) gray (color) scales represent positive and negative toroidal 
fields, respectively; while the continuous and dashed lines rep- 
resent the positive and negative poloidal fields, respectively. For 
this model r=13.6 yr, 5^=4.06 x 10 4 G and B rm> =\46.% G. 
Only toroidal fields larger than 2 x 10 4 (the most external con- 
tours) are shown in panels (a) and (b). This model has started 
with anti-symmetric initial condition (see §5 for details.) 

dynamo which represents a large scale version of the com- 
ponent of the a tensor. It should resemble the buoyant rising and 
twisting of strong magnetic flux tubes, reproducing events of fast 
emergence of them (see the continuous lines in the diagrams of 
Figure 1). We note that another source of poloidal field may ex- 
ist inside the convection zone and at the tachocline, but for our 
purpose we consider only the Babcock-Leighton a effect since it 
is in fact, observed at the surface. 

We solve equation (Q]) for A and with r and 6 coordinates in 
the spatial ranges Q .6R G -Ro and 0-n, respec tively, in a 200 x 
200 grid resolution dGuerrero & Munoz1l2004 see details of the 
numerical model in). 

3. Results 

The velocity field (U) considered in the calculations below 
corresponds to the anal y tical p rofiles of eqs. (4) and (5) of 
iDikpati & Charbonneaul (1 19991) . According to the Babcock- 
Leighton mechanism, the alpha term (aB^,) is concentrated be- 
tween O.95R and R G and at the latitudes where the sunspots 
appear (see the continuous lines in Fig. 1). Since it must pro- 
duce the emergence of magnetic flux tubes, we consider this 
term as being proportional to the toroidal field B,j,{r c ,9) at the 
overshoot interface r c =0.715/? o . For the magnetic diffusion, we 
consider only one gradient of diffusivity located at r c , which sep- 
arates the radiative stable region (with rj^-lQ 9 cm s~ 2 ) from the 
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Fig. 3. The same as in Figs. 2 but considering here the radial 
diamagnetic term. For this model T-13.6 year, 5^=3.9x10 G 
and ,-,„,.= 147.6 G. This model has started with anti-symmetric 
initial condition. 



convective turbulent one (with ?7 cz =10 n cm s~ 2 ) (see the dotted 
line in the upper panel of Fig. 1). The non - dimen sional parame- 
ters as defined in iDikpati & Charbonneaul d!999h and employed 
in the models have the following values: R m -UoR /r] cz -695.5, 
Ca=Q e qR%lric:=\A x 10 5 and C a =aoR Q /ri c: =ll.8, where R m is 
the magnetic Reynolds number, Uq = is the maximum merid- 
ional flow velocity, Q. eq is the angular velocity at the equator, and 
ao is the maximum amplitude of the a effect. In all the cases, the 
tachocline thickness corresponds to 2% of the solar radii and is 
at a radius R = 0.7R Q . 

In Figure 2, the turbulent pumping is not considered. As it has 
been reported in GDPa and GDPb, the most important contribu- 
tion to the toroidal field comes from the latitudinal shear term 
and therefore, the field that is responsible for the observed ac- 
tivity begins to be formed inside the convection zone. On one 
hand, it can be seen that the penetration into the stable layer is 
very weak, on the other hand, the equartorward velocity is faster, 
so that the time that the toroidal field has to amplify to the val- 
ues required by the magnetic flux tube simulations is probably 
short. The upper panel of Fig. 2 shows the butterfly diagram af- 
ter a transient time of 3.3 x 10 4 years. At this time, the toroidal 
field has almost reached a quadrupolar parity (i.e. the toroidal 
field in both hemispheres has the same sign), which is in con- 
tradiction with the Hale's law. In the solar dynamo modelling 
this problem is known as the parity problem and we will discuss 
this subject in more detail in section §5. In the next sections we 
consider the turbulent magnetic pumping as an alternative mech- 
anism of penetration. The latitude of emergence of the toroidal 
field depends on the sta bility criterion for the buoyancy (e.g. see 
iFerriz-Mas et al.|[l994l and their Figs. 1 and 2 for details). 



a) Butterffy diagram 
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Fig. 4. The same as in Fig. 2 for a model with the full turbulent 
pumping terms obtained from magneto-convection simulations. 
For this model T=8.2 yr, B^=4.8 x 10 4 G and B r ,„„=155.92 G. 
This model has started with anti-symmetric initial condition. 

3. 1 . Diamagnetic pumping 

A way to describe the diamagnetic behaviour of a non- 
homogeneous plasma based on a first order smoothing approxi - 
mation (FOSA) was outlined bv lKichatinov & Ruedigerldl992l) . 
If there is an inhomogeneous diffusivity in a fluid, it causes trans- 
port of the magnetic field with an effective velocity. We intro- 
duce this effect in the model by changing the term 7/7- V x B to 
tjt V x B + V777- x B/2 in eq. £[). This new term increases the 
poloidal flux that penetrates beneath the overshoot layer (Figure 
3) and, as a consequence, more toroidal field is produced at all 
latitudes. The toroidal field formed inside the convection zone 
penetrates the stable layer where it is amplified by the radial 
shear at the tachocline. The equartorward velocity of the mag- 
netic flux inside the stable layer is smaller than the velocity right 
above the overshoot interface, so that in the absence of latitudinal 
pumping, the toroidal field would last longer at this layer. The ra- 
dial velocity corresponding to this effect at the overshoot region 
is f/ i /, fl =-V?77-/2^47 cm s _1 when a variation of two orders of 
magnitude is considered for the diffusivity in a thin region of 
0.015fl s . 

3.2. Full pumping 

In a convectivelly unstable rotating plasma, the magnetic field 
is not advected in the vertical direction only. The diamagnetic 
effect may have components in all directions. Besides, another 
pumping effect due to density gradients can develop and in some 
conditions can produc e an upward transport th at can balance 
the diamagnetic effect (Ziegler & Riidiger 2003). Aiming at in- 
vestigating a more general pumping advection, we will con- 
sider in this section the integration of eq. (1) with £ given 
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by eq. (2) and y given by both contributions. The radial and 
latitudinal components of this total y were computed numer- 
ica lly from three-dim e nsiona l magneto-con v ection simulations 
b v lOssendriiver et al.1 d2002l) : iKapvla et all (l2006al) . Similarly 
to IKapvla et al. ( 2006b ), we use the following profiles approx- 
imately fitted from the numerical simulations: 



7e = 7oe 



1 +erf 



r -0.8 



1 - erf 



0.55 
r - 0.98 



(3) 



0.025 



x cos sin 6. 



-70r 



1 +erf 



r- 0.715 



1 - erf 



0.015 
r-0.97^ 



(4) 



0.1 



r- 0.715 , 

:xp |-6-oT^' cos9 + 



where yog and yog define the maximum amplitudes of the pump- 
ing coefficients. The latitudinal pumping is zero at the overshoot 
layer and assumes positive (negative) values at the convection 
zone in the north (south) hemisphere; it is zero at the poles and at 
the equator, with a maximum value of 100 cm s _1 around ~ 15° 
(see the dot-dashed lines in Fig. 1). The radial pumping, y r , is 
negative at the convection zone, indicating a downward transport 
until r=0.7 (i.e., below the overshoot interface) and vanishes be- 
low this value. Its maximum amplitude is 40 cm s _1 (see dashed 
lines in Fig. 1). We do not consider the longitudinal contribution 
of the pumping because it is small compared with the longitudi- 
nal velocity term w<j. 

Figure 4 shows that besides advecting the magnetic fields to- 
wards the stable regions, the pumping terms lead to a distinct 
latitudinal distribution of the toroidal fields when compared with 
the previous results of Figs. 2 and 3. The turbulent and density 
gradient levels present in a convectively unstable layer cause the 
pumping of the magnetic field both down and equartorward, al- 
lowing its amplification within the stable layer and its later emer- 
gence at latitudes very near the equator. This result is important 
for the dynamo modeling because it suggests that the pumping 
can not only solve the problem of the storage of the toroidal 
fields in the stable layer, but it can also help to provide a latitu- 
dinal distribution that is in agreement with the observations. 

3.3. A shallow meridional flow 

As the pumping and the meridional flow are both advective terms 
and in some regions inside the convection zone their radial and 
latitudinal components have the same sign, when the total pump- 
ing is considered the period of the cycle is strongly affected. It 
goes from 13.6 yr in the models of Figs. 2 and 3, to 8.2 yr in the 
model of Fig. 4. One possibility to reproduce the solar period is 
to decrease the value of the diffusivity at the convection zone. 
Another possibility is to decrease the depth of penetration of the 
meridional flow. This is supported by recent helioseismic results 
dMitra-Kraev & Thompsonl (|2007[) ) that suggest that the return 
point of the meridional circulation can be at ~ O.95/? . At lower 
regions, beneath ~ 0.8/? Q , a second weaker convection cell or 
even a null large scale meridional flow can exist. In Figure 5, we 
have decreased the depth of penetration of the flow and found 
that the period increases to the observed value at the same time 




Fig. 5. Meridional flow streamlines and the butterfly diagram for 
a model with the full pumping term, but with a shallow merid- 
ional flow penetration with a depth of only O.8/? , £/o=1300 cm 
s _1 , y«)=90 cm s _1 and y,o=30 cm s _1 . For this model we obtain 
F=10.8 yr, B 0m =4.5x 10 4 G and fi w =154.9 G. This modelhas 
started with anti-symmetric initial condition. 



that the toroidal fields become more concentrated at lower lati- 
tudes. If we further decrease the depth of penetration, the equa- 
torward concentration of the toroidal fields becomes larger and 
the period longer. The reason for that is that no net magnetic flux 
is going poleward at the lower regions of the convective layer 
since it is going all to the equa tor with the pumping vel ocity. 
This result is in agreement with lOssendriiver et alj y()02) who 
suggest than the equartorward motion of the magnetic activity 
could not be the result of a meridional bulk motion, but due to 
the latitudinal pumping of the toroidal mean magnetic field. A 
parametric analysis of the simulations performed under the con- 
ditions above results: 



I 500 < Uo < 3000 cm s" 



T « 181.2 C/ 012 r,fl -'VeT 5 , | 60 < y m < 140 cm s" 1 

[ 20 < y,-o < 120 cm s" 1 



(5) 



This indicates that the pumping terms regulate the period of the 
cycle, leading to a different class of dynamo that is advection- 
dominated not by a deep meridional flow but by turbulent pump- 
ing. Figure 5 shows a butterfly diagram with fiducial values for 
the meridional flow which result a period of 10.8 yr. The agree- 
ment of this diagram with the main features of the solar cycle, 
including the phase lag between the field components, is clear. 

4. Q effect at the near-surface shear layer 

Helioseismology inversions have identified a second radial shear 
layer locat ed below the solar photosphe re in the upper 35 Mm 
of the sun dCorbard & Thompsonll2002l) . It is possible, that the 
solar dynam o is operating in this region, as has been critically 
discussed bv iBrandenburd d2005l) . The more attractive features 
of an aQ. dynamo operating in this region are, among others, 
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Fig. 6. Butterfly diagram for a model with the same parame- 
ters of Fig. 2, but with near-surface shear layer. For this model 
r=15.6 yrB^(r=0.715)=l.l x 10 5 , fi 0mr , t (r=O.98)=1.9 x 10 4 G 
and B,. |rar =131.7 G. This model has started with symmetric initial 
condition. 



the following: (i) the intensity of the magnetic flux tubes in this 
region does not need to be as large as 10 5 G in order to form 
sunspots with the observed magnitudes, but 10 3 G is sufficient; 
and (ii) with near-surface Q effect it is possible to explain the 
coincidence of the angular velocity of the sunspots in the pho- 
tosphere with the r otation velocity at R-0.95R G (see Fig. 2 of 
lBrandenburdl2005l) . as well as t he apparent discon nection be- 
tween the sunspot and its roots dKosovichevI l2002h . The con- 
tribution of a near-surfa ce radial shear has been investigated in 
interface-like dynamos dMason et alj |2002|). in distri buted dy- 
namos with a turbulent a effect dKapyla et al. 2006b), and also 
in advection dominated dynamos dDikpati et al.ll2002l) . The lat- 
ter authors have discarded the radial shear layer since it gener- 
ates butterfly diagrams in which a positive toroidal field gives 
rise to a negative radial field, which is exactly the opposite to the 
observed. In this section, we include the radial shear term in our 
FDBL model in order to explore the contribution of the pump- 
ing to this new configu ration. We use the an alytical expression 
given in eqs. (l)-(3) of lDikpatietai] d2002l) . The near-surface 
shear described by these equations relates a negative shear be- 
low 45° with a pos itive shear above this latitude (see Fig. 1 of 
iDikpati et al.ll2002h R 

With the assumption that the sunspots are formed in the upper 
layers, the Babcock-Leighton poloidal source term, which is 
concentrated in the same region (above r=O.95/? ), does not 
have to be non-local anymore. For the same reason, the values 
of both, the radial and the toroidal fields in the butterfly diagram 
can be taken in the same radial point (r=0.98/? G ). Using the 
same parameters as in the model of Fig. 2, but considering 
a near-surface shear, the results of Figure 6 show two main 
branches in the butterfly diagram. One is migrating poleward (at 
the high latitudes) and one is migrating equatorward (bel ow 45 o. 
This result is expected if the Parker- Yoshimura sign rule dParkerl 
1 1 955b lYoshimural 1975h is considered. We note that the resulting 
parity is quadrupolar but with the correct phase lag between th e 
fields, which is opposite to the obtained in lDikpati et alj d2002l) . 
This difference probably arises from the fact that we are using 
a lower meridional circulation amplitude. Anyway, the polar 
branches are strong enough to generate undesirable sunspots 
close to the poles. The period increases to 15.6 y, which is due 
to the fact that the dominant dynamo action at the surface goes 
in the opposite direction to the meridional flow. 



1 Thi s profile is slightly different from that used by iKapyla et alj 
d2006bh since the latter consider a negative radial shear at all latitudes. 
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Fig. 7. The same as in Fig. 5 for a model with near-surface 
shear action. For this model T=16.3 yr B K rav (r=0.715)=9.7xl0 4 , 
B (r=0.98)=1.9 x 10 4 G and B w =164.4 G. This model has 
started with symmetric initial condition. 

In Figure 7 the same parameters as in Fig. 5 have been used, 
but this time considering the radial shear near the surface. As 
the radial pumping has its maximum amplitude close to the 
poles (see the dashed line in Fig. 1), the toroidal fields created 
there are efficiently pushed down before reaching a significant 
amplitude, so that only the equatorial branches below 45° 
survive. This scenario requires that the pumping be dominant 
over the buoyancy at such latitudes. Also, the phase relation 
of Bf Btf,, obtained in the model of Fig. 7 seems to be the one 
observed, at least at the latitude of activity, however, there is 
some overlapping between one cycle and the next. Results 
which are in better agreement with the observations may be 
achieved if the parameters are finely tuned. 

We note that the introduction of the radial shear close to the sur- 
face when a meridional flow cell penetrating down to r=O.71/? 
is considered, as in the model of Fig. 4, requires an increase of 
th e amplitude o f a. Thi s result is in agreement with that found 
bv lKapvlaetaT]d2006bl) 

5. Brief remarks on the parity problem 

Despite that it has been already explored by several authors, 
the anti-symmetry (dipolar parity) or symmetry (quadrupolar 
parity) of the toroidal magnetic fields across the solar equator 
still constitutes one of the most challenging questions in the 
solar dynamo theory. This is mainly because the resulting 
parity in a model is very sensitive to a huge parameter space. 
The solar-like (antisymmetric) solution could result from 
the effective diffusive coupling of the poloidal field in both 
hemispheres dChatteriee et alj|2004l) . but it may also depend on 
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the position and amplitud e of the a effect (iDikpati & Gilmaril 
l200UlBonanno et al.ll2002l) . or be the result of the imprint of the 
quadrupolar form of the meridional flow on the poloidal mag- 
netic field, as argued by ICharbonneaul (120071) . Small variations 
in the parameter space can switch one solution from a dipolar 
to a quadrupolar one. Although the main goal of this work was 
not to study the parity problem itself, but the contribution of the 
turbulent magnetic pumping, it is interesting to take advantage 
of the full sphere integration in order to see how the pumping 
affects the parity. 

All the simulations presented in the previous sections evolved 
10 7 time steps up to ~10 4 years. All started with antisymetric 
(A) or with symmetric (S) toroidal magnetic field, but we have 
also performed tests with rand om (R) fields. The parity of the so- 
lution is calculated, following IChatterjee et al.l d2004l) . with the 
equation below: 



a) Without pumping 



Pit) 



f-Tn iBl 



J j T l/ 2 (B N (t) - B N ) 2 dt J J T ^ 2 (B s (t) - Bsfdt 



(6) 



where and B$ are the values of the toroidal magnetic field 
at r-Q.7l5R & , and 0=25° and -25°, respectively, and B$ are 
their respective temporal averages over one period. The value of 
P should be between +1 (symmetric) and -1 (anti-symmetric) 
depending on the parity of the fields. The results of our simula- 
tions with regard the parity can be summarized as follows: 

- The models without pumping (e.g., Fig. 2), and those with 
diamagnetic pumping (Fig. 3) result in quadrupolar solu- 
tion. When beginning with a dipolar initial condition they 
take several years before switching to a quadrupolar solu- 
tio n (see Fig. 8a). This result div erges from the on e obtai ned 
by IDikpati & Gilmanl d200ll) or IChatteriee etail d2004l) in 
which the change begins only after around 500 yr. This result 
indicates the strong sensitivity of the parity to the initial pa- 
rameters, in such a way that, for example, the present parity 
observed in the sun could be temporary, at least in the case 
that the turbulent pumping is not relevant for the dynamo. 

- The models with full pumping (e.g., Fig. 5) conserve the ini- 
tial parity if this is symmetric or anti-symmetric (see Fig 8b). 
When initialized with a random field, the system tends first 
to choose the quadrupolar parity, but then it tends to migrate 
to the anti-symmetric (dipolar) parity (dot-dashed line of Fig 
8b), suggesting that the strong quadrupolar imprint due to 
meridional circulation could be whashed out when the full 
turbulent pumping is switched on. 

- Models with full pumping plus near-surface shear layer (e.g. 
Fig. 7) tend to the dipolar parity since the first years of in- 
tegration (see continuous line of Fig. 8c). In these models 
we find that the coupling of the polodial fields in both hemi- 
spheres is more effective. This is probably d ue to the employ- 
ment o f a local a term. See, for example also lChatterjee et all 
(12004 . where a near-surface a effect is combined with a 
buoyancy numerical mechanism. They find a similar cou- 
pling. However this coupling alone is not enough to ensure 
a dipolar parity (see Fig. 6 and the dashed line of Fig. 8c). 
It is also necessary to eliminate the effect of the quadrupo- 
lar shape of the meridional flow upon the poloidal magnetic 
component. This can be done by the action of the pumping 
at the entire convection zone (as indicated by the continuous 
line of Fig. 8c). 
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b) With full turbulent pumping 
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Fig. 8. Parity curves for the three classes of models considered, 
i.e., (a) for models without pumping (as e.g., in Fig. 2); (b) for 
models with full pumping (as e.g., in Fig. 5); and (c) for models 
with near-surface shear (as e.g., in Figs. 6 and 7). In the panels 
(a) and (b), the continuous, dashed and dot-dashed lines corre- 
spond to symmetric, anti-symmetric and random initial condi- 
tions, respectively. In the bottom, (c), panel, the continous line 
is used for the model with turbulent pumping while the dashed 
line is for the model without pumping. 

6. Discussion and conclusions 

We have performed 2D numerical simulations of BLFD solar 
dynamo models including the turbulent pumping. Our first set 
of simulations include a solar rotation profile but without the a 
near-surface radial shear layer. The results show that the pump- 
ing transport effect is, in fact, relevant in solar dynamo mod- 
elling, since it can solve two important problems widely dis- 
cussed in the literature: the storage of the toroidal field at the sta- 
ble layer and its latitudinal distribution. A new class of dynamo 
is proposed in which the meridional flow is important only near 
the surface layer in order to make the Babcock-Leighton mecha- 
nism to operate over the toroidal fields, while in the inner layers, 
the advection is dominated by the pumping velocity. Our results 
support the idea that the equatorward migration of the sunspot 
activity is related to the latitudinal pumping velocity at the over- 
shoot layer and the convection zone. Another attractive feature 
of this model is that a large coherent meridional flow is not any 
more required. 

In a seco nd set of simulations, we have included the shear layer 
found by ICorbard & Thompsonl d2002l) at the upper 35 Mm of 
the sun. The results show the formation of a second shell of 
strong toroidal field just below the photosphere when the full 
pumping is absent. The branches of this field obey the Parker- 
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Yoshimura sign rule for a positive a effect, i.e., they move pole- 
ward at high latitudes and equatorward below 45°. The role of 
the pumping on this kind of models is also interesting since it 
reduces the amplitude of the polar toroidal fields pushing them 
inwards (Fig. 7). It should be noted that these models only work 
fine if a shallow meridional circulation profile is used. When a 
deeper meridional flow going down to the tachocline is consid- 
ered, a strong a effect is required in order to excite the dynamo. 
With regard to the parity problem, our results show that a sim- 
ple cO dynamo with the a effect concentrated near the surface 
leads to a quadrupolar parity, although the switch from dipolar 
to quadrupolar parity takes longer than in previous studies. The 
models with full pumping conserve the initial parity, and when 
the initial condition is random, the system tends to switch to a 
dipolar parity. All the models that combine full pumping with 
near-surface shear prefer dipolar parity solutions too. 
In summary, our results have demonstrated the importance of the 
pumping in the solar dynamo, and suggest that this effect must 
be included in forthcomin g studies, even in those that employ 
multi ple convection cells dBonanno et al.l 120061 : iJouve & Brunl 
120071) . Besides it decrease the influence of the meridional flow 
in two important aspects: on the period of the cycle and on the 
latitudinal distribution of the toroidal fields. On the other hand, 
our results indicate that in the presence of full pumping there 
are two possible solutions to the question on where the dynamo 
operates: it could be either at the convection zone with the mag- 
netic flux tubes emerging from the overshoot layer, or it could be 
at the layers near the surface. Both possibilities have their pros 
and cons, however a kinematic dynamo model only is not suffi- 
cient to conjecture a definitive answer and we will explore this 
in more detail in forthcoming work. 
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